New Type of Magneto-Rotational Instability in Cylindrical Taylor-Couette Flow 
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We study the stability of cylindrical Taylor-Couette flow in the presence of combined axial and 
azimuthal magnetic fields, and show that adding an azimuthal field profoundly alters the previous 
results for purely axial fields. For small magnetic Prandtl numbers Pm, the critical Reynolds number 
Re c for the onset of the magneto-rotational instability becomes independent of Pm, whereas for 
purely axial fields it scales as Pm' 1 . For typical liquid metals, Re c is then reduced by several 
orders of magnitude, enough that this new design should succeed in realizing this instability in the 
laboratory. 
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The magneto-rotational instability (MRI) is one of the 
most important processes in astrophysics, with applica- 
tions ranging from accretion disks to entire galaxies jlLul- 
There is therefore great interest in trying to study this 
instability in the laboratory @,S HUH SUES- These 
experiments are complicated by the extremely small mag- 
netic Prandtl numbers of liquid metals, which necessitate 
very large rotation rates. Here we report on a new type 
of magneto-rotational instability, that operates even at 
infinitesimal Prandtl number, and correspondingly at far 
lower rotation rates. This new instability should be much 
easier to obtain experimentally, and may also have con- 
siderable astrophysical implications. 

The magneto-rotational instability is a mechanism for 
transporting angular momentum. Consider, for example, 
an accretion disk around a newly forming star. In order 
for material to fall inward, it must transfer its angular 
momentum to material further out. The difficulty is how 
to accomplish this. In particular, a Keplerian angular 
velocity profile, for which £1 ~ r -3/2^ j g i cnown ^ Q b e n y_ 
drodynamically stable (at least linearly), by the familiar 
Rayleigh criterion. Purely viscous coupling in a laminar 
flow is many orders of magnitude too small though; if 
this were the angular momentum transport mechanism, 
stars would take so long to form that none would yet 
exist today, some 14 billion years after the Big Bang. 

This conundrum was solved by Balbus & Hawley 0, 
who showed that such a Keplerian profile may be hydro- 
dynamically stable, but is nevertheless magnetohydro- 
dynamically unstable. The addition of a magnetic field 
opens up a new way of coupling fluid parcels, namely 
via the magnetic tension in the field lines, and thereby 
allows for a new instability, the magneto-rotational in- 
stability, that has no analog in the purely hydrodynamic 
problem. Coupling the fluid in this way is then so much 
more efficient at transporting angular momentum that 
this mechanism can indeed yield accretion rates more in 
line with those observed. Something as basic as the time 
it takes for a star to form is thus magnetically controlled. 



Subsequent to this work, it was soon realized that this 
instability had actually been discovered decades earlier, 
by Velikhov wn0 na cl not considered it in an as- 
trophysical context though. Instead, he had viewed it 
simply as the magnetohydrodynamic extension of the so- 
called Taylor-Couette problem, consisting of the flow be- 
tween differentially rotating cylinders. See for example 
[lij l for a review of this problem, one of the oldest in fluid 
dynamics. 

Once the connection is made to the classical Taylor- 
Couette problem, one immediately recognizes that this 
should be a way to study this instability in the lab: just 
take the fluid to be a liquid metal, and apply a mag- 
netic field along the cylinders. This simple design is the 
basis of most of the MRI experiments proposed to date 
dH HIE [3,0 (Sisan et al. @ claim to have obtained 
the magneto-rotational instability already, in spherical 
rather than cylindrical geometry. However, there are 
enough other instabilities that can arise in this config- 
uration, e.g. [12} i that their results are also open to 
other interpretations. Indeed, the basic state from which 
their instability arises is already fully turbulent, indicat- 
ing that whatever their instability may be, it is certainly 
not the first instability to set in.) 

As simple as it sounds, there is unfortunately also 
one considerable difficulty associated with these exper- 
iments, namely that the rotation rates of the inner and 
outer cylinders must be enormous. The problem is that 
the relevant parameter is not so much the hydrodynamic 
Reynolds number Re = fij?"f jv, but rather the magnetic 
Reynolds number Rm = flirf /rj, where v is the viscosity 
and rj the magnetic diffusivity. In its traditional form, the 
magneto-rotational instability sets in when Rm = O(10). 
Re is then given by Rm/Pm, where Pm = vjr\ is the 
magnetic Prandtl number, a material property of the 
fluid. Typical values are ~1CT 5 for liquid sodium, ~10~ 6 
for gallium, and ~10~ 7 for mercury. Re must therefore 
be at least 10 6 , and possibly larger still, depending on 
the particular liquid metal one intends to use. 
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Such large values can be reached in the lab; taking the 
inner cylinder radius rj = 10 cm, say, one finds that ro- 
tation rates of around 10 Hz are required - large, but 
achievable. The difficulty lies elsewhere, namely in the 
ends that would necessarily be present in any real ex- 
periment. Due to the Taylor-Proudman theorem, stat- 
ing that in rapidly rotating systems the flow will tend 
to align itself along the axis of rotation, end-effects be- 
come increasingly important, until at Re — O(10 6 ) the 
flow is controlled almost entirely by the end-plates 
That is, while it is possible to achieve such large rota- 
tion rates, the flow will look nothing like the idealized 
infinite-cylinder flow on which all of the theoretical anal- 
ysis is based. 

A new approach is therefore needed if this instability 
is to be studied in the laboratory. We propose the fol- 
lowing: instead of imposing only a uniform axial field 
B z = Bo, additionally also impose an azimuthal field 
£?0 = f3Bo(ri/r), where (3 is a nondimcnsional measure 
of the relative magnitudes of B^ and B z , and the dimen- 
sional quantity Bq will be incorporated into the Hart- 
mann number below. Such a B^ field can be generated 
just as easily as a B z field, by running a current-carrying 
wire down the axis of the inner cylinder (suitably insu- 
lated from the fluid, of course). Note also that with nei- 
ther B z nor B^ maintained by currents within the fluid 
itself, the possibility of magnetic instabilities is excluded 
a priori. The only source of energy to drive an instabil- 
ity is the imposed differential rotation; the magnetic field 
merely acts as a catalyst. 

Given this basic state consisting of these externally im- 
posed magnetic fields, as well as the differential rotation 
profile f2(r) imposed by the rotation rates fii and fl Q of 
the two cylinders, we linearize the governing equations 
about it, and look for axisymmetric disturbances. These 
are known to be preferred over non-axisymmetric ones 
for the classical (3 — MM, and are therefore the appro- 
priate starting point here as well. The perturbation flow 
and field may then be expressed as 

u = ve^ + V x (tpettj,), b = be^ + V x (ae^). 

Taking the z and t dependence to be exp(ikz + "ft), the 
perturbation equations become 

Rejv = D 2 v + Reik r" 1 (r 2 0)' i/j + Ha 2 ik b, 

Re 7 D 2 i\) = L> 4 V - Re2ikflv + Ha 2 ik(D 2 a + 2f3r~ 2 b), 

Pm Re 76= D 2 b — Pm Re ik Vl'r a + ikv — 2ik (3r~ 2 i\), 

Pm Re 7 a ~ D 2 a + ikip, 
where D 2 = V 2 — 1/r 2 , and the primes denote d/dr. 



Length has been scaled by fj, time by f2^~ , by fit, 
u by rj/ri, and Bo and b by Bq. The various nondimcn- 
sional parameters are then: (a) the magnetic Prandtl 
number Pm = v/rj already mentioned above, (b) the ra- 
tio jl — fl /£li (this enters into the details of f2(r) = 
ci +c 2 /r 2 ) and the Reynolds number Re — Qir 2 /i>, mea- 
suring the relative and absolute rotation rates of the two 
cylinders, (c) the parameter (3 and the Hartmann number 
Ha = B ri / y/pjlrji* (p is the fluid's density, \i the mag- 
netic permeability), measuring the relative and absolute 
magnitudes of the imposed magnetic fields, and (d) the 
radius ratio ri/r a , which we fixed at 1/2. 

The radial structure of v, ?p, b and a was expanded 
in terms of Chebyshev polynomials, typically up to N — 
40 — 80. These equations and associated boundary con- 
ditions (no-slip for u, insulating for b) then reduce to a 
large (4iV by 4iV) matrix eigenvalue problem, with the 
eigenvalue being the growth or decay rate 7 of the given 
mode. Note also that this numerical implementation is 
very different from that of 4], in which the individual 
components of u and b were used, and discretized in r 
by finite differencing. Both codes yielded identical re- 
sults though in every instance where we benchmarked 
one against the other. 

The following sequence of calculations was then carried 
out. First, we fixed Pm, /3 and jt, and scanned through a 
range of values of Re, Ha and k, in each case determining 
whether the given modes grow or decay. We thus found 
the smallest value of Re, and the corresponding Ha and 
k, that yields a marginally stable solution, one having 
Re(7) = (which typically involved solving the basic 
eigenvalue problem for several thousand combinations of 
Re, Ha and k). By repeating this entire procedure for 
different values of Pm, (3 and ji, we obtained the results 
in Fig. 1. 

The dotted curve to the left of jl — 0.25 is a purely hy- 
drodynamic instability, namely the onset of Taylor vor- 
tices, e.g. 0. As we approach jl — 0.25 though, we note 
that Re c — > 00. This critical value ji = 0.25 ((f*i/r Q ) 2 , in 
general) is precisely the so-called Rayleigh line, beyond 
which the flow is hydrodynamically stable, because the 
angular momentum increases outward. 

We are more interested therefore in the behavior to the 
right of ji = 0.25, where we anticipate that the inclusion 
of magnetic effects will yield the MRI. We begin with the 
two curves (3 — 0, the MRI as it has been considered to 
date [Tjj. We see how Re c rises very steeply from the 
previous nonmagnetic results to the left of jl = 0.25, and 
then scales as Pm^ 1 , exactly as described above, and 
in perfect agreement with 0. And again, because Pm 
is so small, these values end up being too large for the 
experiment to succeed Q. 

However, turning next to the results for non-zero f3, 
we see that Re c is dramatically reduced, over a range 
of jl extending increasingly far beyond the Rayleigh line. 
For example, if we focus on how far we can go before 
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FIG. 1: The critical Reynolds number Re c as a function of fx. 
The dotted curve to the left of fi — 0.25 is a purely nonmag- 
netic instability. The curves to the right of fi = 0.25 are the 
magneto-rotational instability. Solid curves are Pm = 10~ 6 , 
dashed curves Pm — 10 -5 . The numbers beside pairs of 
curves indicate fi. Finally, the dot on the f3 — 4 curves corre- 
sponds to the solution shown in Fig. 2. 

Re c = 10 4 , say, we obtain ft = 0.253, 0.264, 0.292 and 
0.308 for f3 = 1, 2, 4 and 8, respectively. Furthermore, 
within these fi ranges Re c is independent of Pm, very 
different from the previous Pm -1 scaling, indeed, within 
these ranges one can set Pm = and still obtain exactly 
the same solutions. This Pm = limit was considered 
before by Chandrasekhar but for axial fields only, 
in which case there are no instabilities to the right of the 
Rayleigh line. 

The explanation for this very different behavior for 
non-zero (3 lies in the coupling between the azimuthal 
field b and the meridional circulation ip. If (3 = these 
quantities are not directly coupled at all, only indirectly 
through v and a. For non-zero f3 each acts directly on 
the other. Of these two new terms in the equations, the 
more important one turns out to be the effect of if on 
b. Physically, this corresponds to the meridional circula- 
tion if advecting the original and thereby generating 
a contribution to b. ft is this new way of maintaining b 
that allows this instability to proceed even in the Pm = 
limit, where the classical (3 — MRf fails. 

Figure 2 shows an example of these new solutions. We 
note that however dramatically Re c may have been re- 
duced, the spatial structure is much the same as for 
(3 = 0, consisting of Taylor vortices elongated slightly 
in the z-direction. There is actually one subtle differ- 
ence, namely that the up/down symmetry in z has been 
broken. This is due to the handedness of the imposed 
field, which distinguishes between ±z in a way that a 
purely axial field does not. As a result of this symmetry- 
breaking, these new modes are also no longer station- 
ary, but instead have fm(7) ^ 0, corresponding to a 
drift in z, at the rate Im("/)/k. This particular solu- 
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FIG. 2: The spatial structure of the eigenmode for Pm — 0, 
fi = 4 and jj, = 0.27, for which Re c = 1521 (see Fig. 1), 
Ha — 16.3 and k = 2.33. From left to right, contours of v, ip, 
b and a, with the contour interval for b and a 1/3 that of v and 
ip. Note the lack of any particular phase relationship between 
the various quantities, due to the ±z symmetry-breaking. 



tion has Im(7) = 0.153 and k = 2.33, for a drift rate of 
0.066 (rA)- 

Turning next to the Hartmann number Ha = 16.3, this 
translates into a ~12 G field, taking the fluid to be liquid 
sodium, and rj = 10 cm. An axial field of that magni- 
tude is certainly easily achievable in the lab. For the 
azimuthal field we then want 500 G cm/r, corresponding 
to a current of 2500 A in this wire running down the cen- 
tral axis, which is again achievable, if perhaps not quite 
so easily. 

To summarize then, we have shown that the magneto- 
rotational instability is radically altered if one imposes 
both axial and azimuthal magnetic fields, becoming inde- 
pendent of the magnetic Prandtl number in the Pm — > 
limit - a result that is all the more remarkable as neither 
purely axial nor purely azimuthal fields yield anything 
like it. With this new scaling, the critical Reynolds num- 
bers are reduced by several orders of magnitude, enough 
that these instabilities could perhaps be attainable in the 
lab without being disrupted by end-effects (although of 
course some end-effects will always be present, partic- 
ularly with this slow drift in z). Further computational 
work includes the nonlinear equilibration of these modes, 
as well as the possibility of non-axisymmetric instabili- 
ties. These and other issues are currently being explored. 

Finally, returning briefly to the original astrophysical 
motivation, we note that virtually all astrophysical bod- 
ies do in fact have both axial and azimuthal magnetic 
fields. These magnetic fields are typically not externally 
imposed though, but rather generated by electric cur- 
rents flowing within the system itself. Sclf-consistently 
solving for both the large-scale fields as well as the small- 
scale instabilities is then far more complicated than our 
analysis here, but the basic ingredients are certainly there 
for this new type of magneto-rotational instability to play 
a role. 
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